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Abstract 

Clan concept has been introduced in multiparticle dynamics in or- 
der to interpret the wide occurrence of negative binomial (NB) regu- 
larity in n-charged particle multiplicity distributions (MDs) in various 
high energy collisions. 

The centrality of clan concept led to the attempt to justify its occur- 
rence within a statistical model of clan formation and evolution. In 
this framework all thermodynamical potentials have been explicitely 
calculated in terms of NB parameters. Interestingly it was found that 
NB parameter k corresponds to the one particle canonical partition 
function. 

The goal of this paper is to explore a possible temperature T and vol- 
ume V dependence of parameter k in various classes of events in high 
energy hadron-hadron collisions. 

It is shown that the existence of a phase transition at parton level 
from the ideal clan gas associated to the semihard component with 
k > 1 to the ideal clan gas of the hard component with k < 1 implies 
a discontinuity in the average number of particles at hadron level. 
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1 Introduction. A recollection of statistical 
and quantum theoretical ideas in multipar- 
ticle production 

One of the most interesting experimental facts in the thierties of the past 
century had been the discovery in cosmic ray physics of multiparticle pro- 
duction from a primary hadron of given energy. 

E. Fermi [1] in order to explain the new phenomenon proposed a thermody- 
namical approach to it. In the same years L. Landau [2] focused his atten- 
tion on the hydrodynamical aspects of the production process, whereas W. 
Heisenberg [3] was impressed by its non-linearity. Accordingly, Heisenberg 
took a different point of view from Fermi and Landau and argued that the 
natural framework for the description of multiparticle production should be 
found not in statistical mechanics but in a non-linear field theory of a "new 
nuclear force" . 

In between these two quite extreme approaches, the first one inspired by sta- 
tistical mechanics and the second one by quantum field theory, are contained 
all the later developments of multiparticle production in high energy physics. 
Today we know indeed that the non-linear quantum field theory foreseen by 
Heisenberg is Quantum Chromo Dynamics (QCD), the quantum field theory 
of the strong interactions (the "new nuclear force"...). 

At the same time the discovery in the sixties of the past century of unsus- 
pected regularities in final particle multiplicity distributions originated by 
primary hadrons with different initial energies when plotted vs. n in cosmic 
ray physics by P.K. Mac Keown and A.W. Wolfendale [I] pointed out the 
statistical aspects of collective variables properties like normalized cumulants 
(K n ) and factorial moments (F n ) of the final particle multiplicity distribu- 
tion describing the production process. It was found in fact that n final pion 
multiplicity distributions, P n , when plotted vs. n turned out to be quite well 
fitted by a NB (Pascal) MD with the average multiplicity, n, and k — D 2_ n 
(D is the dispersion of the MD) as characteristic parameters. In addition, 
n is becoming larger and parameter k smaller as the energy of the primary 
hadron increases. 

One of the main results along this line of search - it should be mentioned 
- had been in view of the decreasing of parameter k, as the primary energy 
increased, clearly KNO scaling violation. 

Once again cosmic ray physics findings anticipated future discovery in the 
accelerator region. 

It has to be noticed in fact that the study of collective variables properties 
in various classes of high energy collisions at various cm. energies and in dif- 
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ferent (pseudo)rapidity intervals in the accelerator region performed in the 
seventies and in the eighties by various Collaborations revealed in pp, pp, \ip 
and amazingly in e + e~ annihilation the same general trends and the same 
regularities discovered in cosmic ray physics. 

In order to reproduce these experimental facts at final hadron level starting 
from perturbative QCD at parton level it is well known that one has to rely 
in a region where QCD has no predictions on Monte Carlo event generators 
based on DGLAP evolution equations with an ad hoc hadronization prescrip- 
tion. This last line of search has to be considered not fully satisfactory in 
view of the quite large number of parameters involved in the event generator 
models and suggests to dig more carefully in the statistical mechanics aspects 
of the approximate universal regularities discovered in P n vs. n plots with 
a limited number of well identified physical parameters. The goal in this 
framework is to understand within the universality of the NB (Pascal) reg- 
ularity the statistical counterpart of the universality of QCD in all explored 
high energy collisions. 

Accordingly, it is instructive to start with a critical reading of the main sta- 
tistical mechanics literature on multiparticle production in addition to the 
pioneering work of Fermi and Landau. 

Among them Feynman fluid analogy should be mentioned. As discussed in 
[5] here the cross-section for the production of n final particles plays the 
role of the partition function in the canonical ensemble and it is an integral 
over phase space of the square of a matrix element, which corresponds to 
the Gibbs distribution e~ H l k B T with Hamiltonian H, Boltzmann factor ks 
and temperature T; volume V is identified by Feynman with the extension 
of phase space and fugacity z with the dummy variable u appearing in the 
definition of the n-particle multiplicity p{n) generating function G{u) 

G{u) = ^u n p{n). 

n 

It is to be remarked that this idetification is again rather unsatisfactory from 
our point of view because one has to satisfy at the same time the definitions 
of the aveage number of particles from the grand canonical ensemble and 
from the definition of the generating function. These definitions can be sat- 
isfied at the same time only in the limit of zero chemical potential. 
Another approach was proposed by Scalapino and Sugar: they defined the 
probability amplitude to produce a particle at rapidity y, denoted by H(y), 
as a random field variable, then they introduced a functional F[U] which 
played a role analogous to the free energy for a system in thermal equilib- 
rium. One can then obtain the n-particle inclusive distribution by averaging 
the product of the squares of the amplitude, II 2 (yi) ■ ■ ■ U 2 (y n ), with a weight 
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given by e~ F ^. Lacking the knowledge necessary to calculate F[U] from the 
underlying dynamics, the authors parametrized it in retaining the first three 
terms in a series expansion, then solved the problem in few particular cases. 
Remarkably, to leading order in the size of the allowed full rapidity range, 
they obtain a generating function of an infinite divisible distribution (IDD), 
a class of distributions to whom the NB (Pascal) MD is belonging. 
More recent results concerning KNO scaling violation and phase transitions 
can be founf in the literature (see Ref.s quoted in [5]). 

Stimulated by these findings some of the present Authors in collaboration 
with S. Lupia proposed a simplified approach to multiparticle production in 
high energy collisions based on the NB (Pascal) MD and valid for any chem- 
ical potential. The interpretation of the NB regularity was given in terms of 
clan structure analysis (see Ref. [S] and Section 2 of the present paper). 
The plan of the paper is the following. 

After introductory remarks discussed in Section 1, in Section 2 the interpre- 
tation of the NB (Pascal) MD regularity in multiparticle production in terms 
of clan structure analysis and its more recent applications to the components 
(or classes of events) in high energy collisions are reviewed. 
In Section 3 a summary of previous results on clan thermodynamics is given. 
In Section 4 the explicit dependence on temperature T and volume V - not 
shown in [5] - of the one particle canonical partition function k, i.e. of the 
NB parameter k in the various components (soft, semihard and hard) of pp 
collisions, is determined. 

In Section 5 it is shown that the existence of a phase transition at parton 
level from the ideal clan gas associated in pp to the semihard component with 
k > 1 to the ideal clan gas associated to the hard component with k < 1 im- 
plies a discontinuity in the average number of particles of the two mentioned 
components at final hadron level. 

In the Appendix a theorem on the zeroes of a degree M polynomial used in 
Section 5 is proved. 

2 Clan concept and Negative Binomial (NB) 
regularity in high energy multiparticle pro- 
duction 

In order to interpret the wide occurrence of negative binomial n charged par- 
ticle multiplicity distribution, P„ , regularity in various classes of high en- 
ergy collisions in the eighties it has been proposed an alternative parametriza- 
tion to standard NB parameters, i.e. the average charged multiplicity n and 
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parameter k of P^ B ^ = fc ( fc+1 )--( fc+ra ~ 1 ) ^_B_^ n ^_A_^ fe The proposed new pa- 
rameters were the average number of groups of particles of common ancestor 
(called in high energy multiparticle phenomenology average number of clans, 
N), and the average number of particles per clan, n c . It has been also shown 
that the new parameters are linked to the standard ones by the following not 
trivial relations 

N = k In ( 1 + - ) and n c = n/N. (1) 



k, 

Notice that n c is a scaling function of n/k. 

It turns out that clans in various classes of events in a collision are indepen- 
dently produced whereas particles within each clan are distributed according 
to a logarithmic multiplicity distribution and all multiplicity correlations are 
exhausted within each individual clan. These properties can be expressed 
by using the corresponding generating functions of the NBMD in the two 
mentioned sets of parameters (n, k) and (N,n c ). In fact 



[NB] x n 



n=0 



k + nil — x) 



(2) 



and 



G^ NB \N,n c ;x) = ex P [N(G ln (n c ;x) - 1)] (3) 
with Gi n (n c ; x) = and z = 

In addition it should be remembered that 1/k is the ratio of the probability 
to have two particles within the same clan over the probability to have in a 
statistical framework the two particles in two distinct clans, i.e. 1/k is an ag- 
gregation parameter which is linked to the two-particle correlation function 
in rapidity variables y 1 and y 2 , C 2 (yi,y 2 ), by \ = J C 2 (yi,y 2 ) dy 1 dy 2 . 
The introduction of clan concept within a single NBMD allowed to single 
out the main differences in multiparticle production in e + e~ annihilation, 
deep inelastic scattering and hadron-hadron collisions. Experts started to 
talk of NB universality in multiparticle production. The regularity at final 
hadron level was supported also by QCD in view of the fact that NBMD 
occurs at the parton level in the leading log approximation with a fixed cut- 
off regularization prescription in order to describe quark- and gluon-jets as 
Markov branching processes. Notice that at parton level 1/k is the ratio of 
the g — > g + g vertex over the q — > q + g vertex, and in the (N, n c ) NB 
parametrization N corresponds to the average number of bremsstrahlung 
gluon-jets, as shown in Ref. [B]. 

After initial enthusiasm which led to a huge list of sound phenomenological 
papers in multiparticle production in high energy collisions, it was found by 
UA5 collaboration at CERN pp collider that the regularity was violated at 
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top cm. energy (900 GeV) and by DELPHI collaboration in e + e~ annihila- 
tion at LEP cm. energy: a shoulder structure occurred in the total n charged 
particle multiplicity distribution, P n , when plotted versus n, which a single 
Pn was unable to reproduce. Interestingly the regularity violated in the 
total sample of events in both reactions was restored at a deeper level of 
investigation, i.e. in terms of its components (or substructures or classes of 
events) general properties: it was found indeed that the regularity initially 
discovered in the low cm. energy domain and then violated in the total sam- 
ple of events at higher cm. energy, continued to be valid in the new energy 
region when applied to the various i-c\\ asses of events (i = 1,2,3...) with 
characteristic ki or iVj, n Cti parameters. Furthermore, it was found that 
the weighted superposition of the various i-classes of events observed in pp 
collisions and in e + e~ annihilation, whose n, particles multiplicity distribu- 
tion was described by a P^ B ^ with characteristic and in general different NB 
parameters for each 2-class of events, allowed to interpret not only shoulder 
structure in the total sample of events when plotted versus n, but also exper- 
imentally observed H q oscillations when plotted versus q (H q is the ratio of 
g-factorial cumulants, K q , and g-factorial moments, F g , of the corresponding 
total multiplicity distribution) as well as forward-backward n-particle multi- 
plicity correlation strenght behaviour in the total sample of events. 
All these facts stimulated the attempt to justify the introduction of clan 
concept and related NB regularity in the various classes of events occurring 
in different high energy collisions in a thermodynamical framework. Clan 
thermodynamics was then proposed [5]. 

3 A summary of previous results on clan ther- 
modynamics 

In paper [5] starting from the relationship between the grand canonical parti- 
tion function Q(z, T, V) of a system with fugacity z, volume V, temperature 
T and the n-particle (parton) canonical one, Z n , for the Negative Binomial 
Multiplicity Distribution (NBMD) with the average multiplicity n and k as 
characteristic NB parameters, it was found that 



oo 



( 



k 



) 



—k 



Q{z, T,V) = z n Z n (T, V) 



n + k 



(4) 



where 



k R T _ 



n 



(5) 



z = e"B 



n + k 
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with fx the chemical potential, Z n = fc ( fc + 1 )--^( A: + w x ) anc [ ^he probability of 
finding n particles (partons) p(n) in the system is given by 

l^n=0 Z Z "A J ' V ) 

with p(0) = {-^£) k and Z = 1. Notice that particle (parton) fugacity z 
satisfies the condition < z < 1, i.e. our particle (parton) system behaves as 
a real boson gas. It should be pointed out that in our approach, as in those 
mentioned in the Introduction, we are always working in the framework of 
equilibrium thermodynamics. 

It turns out that the canonical partition function for a system of one par- 
ticle Zi corresponds exactly to parameter k, i.e. to an unknown function of 
temperature T and volume V 

k = k(T, V). 

Since the thermodynamical potentials calculated in Appendix A of Ref. [S] 
have been shown to depend on the partial derivative of k with respect to 
temperature T, the main goal of the present paper is the determination of 
the explicit dependence of k on temperature T and volume V. 
This search is interesting also in view of the fact discussed in Ref . [7] that NB 
parameter k controls the geometry of the zeroes distribution in the complex 
z-fugacity plane of the M-order algebraic equation of the M-truncated grand 
canonical partition function 

M 

Q M (z, T,V) = Y, znz n(T, V) = 0. (7) 

n=0 

In fact, it has been shown in Ref. [7j that in the rescaled variable u = zx for 
k larger than one all zeroes stay in the complex rescaled w-plane inside the 
circle of unit radius \u\ = 1, whereas for k smaller than one the zeroes stay 
outside the circle \u\ — 1 and on the circle itself for k — 1. 
Accordingly, 

In Q(z, T, V) = N (8) 
or in terms of the Boltzmann equation of state 

PV = k B T In Q = k B TN (9) 

with P the pressure of the ideal gas of clans and k B the Boltzmann constant. 
Volume V is here the space volume occupied by N clans of the ideal clan gas 
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satisfying Boltzmann equation of state (9). 

Notice that here we consider clans and not particles (or partons) as in Eq.(5). 
Accordingly, our system can be interpreted as a real boson gas with fugacity 
z given by Eq. (5) when we are dealing with particles (partons) or as an ideal 
Boltzmann gas when we are dealing with clans. 

These findings have an intriguing application to the eventual new class of 
hard events foreseen at LHC when N — > 1 [5]. In this class of events k is less 
than one, whereas it is larger than one for the soft and semihard classes of 
events. The point to be cleared up concerns indeed the possibility that the 
clan gas in going from k larger than one (i.e. from the semihard clan gas) to 
k smaller than one (i.e. to the hard clan gas) undergoes a phase transition 
and the nature of the phase transition itself. 

Notice that results in the hadronic (h) sector of the collision can be extended 
to the partonic (p) one via the generalized local parton(p)-hadron(/i) duality 
(GLPHD), i.e. by assuming on phenomenological grounds that k p = kh and 
n h = an p with a ~ 2 (see Ref. [9] and related bibliography). 
The existence of the first two classes of events in e + e~ annihilation (class I: 
2-jet events; class II: 3- or more-jet events) as well as in pp collisions (class 
I: events without minijets; class II: events with minijets) was motivated by 
the discovery in both reactions of three relevant facts extensively discussed 
in Ref.s [10] and pj] and mentioned in the introduction of the present paper: 

a. the shoulder structure in n-particle multiplicity distribution of the total 
sample of events P n when plotted vs. n; 

b. H q vs. q oscillations, with H q the ratio of g-factorial cumulants, K q , and 
the (/-factorial moments, F q , of the corresponding MD; 

c. the lack of forward (F) to backward (B) n-particle multiplicity correlation 
in e + e~ annihilation and its strong enhancement in pp collisions. 

It is to be pointed out that all these collective variable properties were ex- 
plained in terms of the weighted superposition of different classes of events, 
each described by a NB (Pascal) MD with characteristic NB parameters. 
The explanation of point c. is very probably a consequence of a huge colour 
exchange mechanism occurring at parton level in pp collisions and of a very 
weak colour exchange in e + e~ annihilation. Clan structure analysis in terms 
of the new parameters N and n c revealed in this framework all its remark- 
able implications in the dynamics of a collision: clans turned out to be more 
numerous in e + e~ annihilation but less populated, whereas just the opposite 
occurred in pp collisions; here clans are less numerous than in e + e~ annihi- 
lation, but more populated. An intermediate situation was shown to happen 
in deep inelastic scattering, i.e. the average number of clans is less numerous 
than in e + e~ annihilation but the average number of particles (partons) per 
clan larger. 
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In conclusion, the dynamics of the collision in terms of N and n c was un- 
derstood as a two step process: clans were firstly independently produced 
(each clan contains at least one particle or parton, its ancestor), each clan 
then decay accordingly to parton or hadron cascading mechanism with all 
particles (partons) correlations exhausted within each clan. 
All these facts led to the conjecture that introduced clan concept is not only 
the way to a new useful and pictorial parametrization to be added to the 
standard one in terms of k and n, but it points in the direction of clans as 
true physical objects to be looked for in future experimental work. In addi- 
tion, the reduction with the increase of the cm. energy of N serni h ard observed 
in pp collisions [12j arises interesting questions on the properties of the even- 
tual third class of quite hard events, characterized by N hard — > 1, k hard <C 1 
and n hard very large with respect to n soft and n semihard . 
By assuming that the third class of events is blowing up already at LHC 
cm. energy (14 TeV), we distinguish at such energy in the collision at least 
three classes of events whose weighted superposition is describing collective 
variable behaviours (i.e. of P n , K n , F n , FB multiplicity correlations) in the 
full sample of events. 

The three expected classes of events in pp are therefore (see Ref.[13j): 

I. a soft one (events with no minijets) due to one parton-one parton scatter- 
ing; 

II. a semihard one (events with minijets) due to two partons-two partons 
scattering; 

III. a very hard one originated by three partons-three partons scattering, 
with N hard — > 1, 

khard 1 and Tlhard ^ ^semihard > f^soft- 

Notice that according to Ref . [13] temperatures as well as the parton densities 
of the i-classes are expected to increase with the hardness of the collision. 
The reduction of the average number of clans as the temperature increases 
(and the parton collisions become more numerous) suggests indeed the start- 
ing of a higher parton density process. 

In the present paper we want to explore the onset of the hard component 
at parton level from the semihard one in pp, assuming it occurs at LHC in 
the framework of clan thermodynamics discussed in [5j, and in particular to 
determine the T and V dependence of NB parameter k. 
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4 Temperature and volume dependence of the 
one particle canonical partition function Z\ = 

In order to calculate analitically the temperature dependence of NB param- 
eter k, we decided to use consistency relation between the entropy expressed 
in terms of the grand canonical partition function Q 

S = k B lnQ + k B T In Qj , (10) 

(/i is here the chemical potential, which is linked to fugacity z as already 
pointed out by z = e k B T ) and the entropy given in the Appendix of Ref . [5] 

where 



and 



F = nk B T\n ( ^— ] - k B T k In (l + ^) 

\n + k J V k / 



N = \nQ = -Hn(l - z) 



It should be noticed that partial derivative with respect to the tempera- 
ture T at constant fi and V can be written in terms of fugacity z as follows 

By asking that Eq.s(10) and (11) coincide, it follows that 

i.e. 
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By applying finally the operator described by Eq.(12) to Eq.(13) and recalling 



that z = we find 



—z In z 



dlnQ 

dz 



+ nlnz = —zlnz 



— z In z 



dz 



(n + k)- ln(l - z) 



dk\ 
Yz) 



+ n In z 



dk 

dz 



+ n In z 



' dk 

zlnzlnil -z) I — - 
1 dz 



N (dk 
J [df 



and 



zlnz 



dk 

dz 



T 



-\n(l-z)T 



dk\ 
df) 



dk 
df 



(15) 



(16) 



Being quantum effects in a boson system with fixed volume V (i.e. constant n) 
a low temperature phenomenon with the fugacity of the system z = =^ — > 1 , 
i.e. k -> 0, one has (^) y > 0. It follows from Eq.(16) that also (||) y > 0. 
In addition k is a canonical variable as it has been shown in [H] and it cannot 
depend on a grand canonical variable like the chemical potential \x. 
Now, willing to extract k variable from (||) y in view of the canonical nature 
of k, and to suppress the zlnz term in Eq.(16), we are led from 



fdk 



dz 



(17) 



nV 



with constant a, to assume that k (lnz) a is a function depending on volume 
V and chemical potential \x only. 
Accordingly, from (17) one has 



dk\ 
dz) 



'dk 



QnzT 



dz 



k (lnz) a 



'd 



(Inzy 

dz 



ak 



zlnz 



> 



with a > 0. Finally from (16) we find that T dependence of k is given by 

(19) 



t(^) =ak. 



\dT 
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Next step is to determine V dependence of k at fixed T. 
By considering now the Helmoltz free energy F 

F = fin — PV 

i.e. in case of the NB 

F = k B T n\nz- k B TN, (20) 

it follows 

(¥) =- p 

which, being in this case PV — k B TN, implies 
f)F\ 1 1 

W ) _ T = -yksTN = —k B Tk ln(l - z). (21) 
Since now from Eq.(20) one has 

n,T 



Z \ dV Jn,T V dV 



_ 2 \( 1 \ 2 f dk\ .,„,,„ N ( dk \ 



~ k ° Tn - 2 {—k [w T +kBTIn(1 - z) {w - 



-k B T(n + k) 



(dz_ 



n,T 



it follows that 

fdF\ , _ f dk\ , ^, , , / dk\ , , s (dz 

W)^~ kBTz w)r +kBTH1 ~ z) W) T ~ kBT(n+k) W 



- - k ° Tz (w) T + ™ " 2 » (§) T " *" r <" + fc » (ll?) 

=-^K§) T +feTin(i - 2) (S, + ^ 2 ®, = 

= k B THl-z) — ^ 



n, 



\dVJ T 



(22) 
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The request that Eq.s(21) and (22) coincide, leads to 

which together with Eq.(19) gives 

k = c T a V 



(23) 



(24) 



with c a constant independent on temperature T and volume V, but depend- 
ing on a only, i.e. a parameter which will be identified - as will be shown 
in the following - with the ideal clan gas degrees of freedom. It should 
be remarked that above results include the ultra-relativistic well of 

course the non-relativistic one, whereas the relativistic case is excluded in 
the present analytical approach and postponed to future work. 
Notice that Eq.(24) is compatible with large T values and N — > 1 (i.e. an 
harder collision process with k — > 0) when the volume is becoming smaller 
(V — > 0) in a faster way than T~ a as the cm. energy y/s becomes larger. 
Equation (24) allows of course to determine all thermodynamical potentials 
given in [5] in terms of c, T and V. It follows: 
Entropy S 



S 



N, 

T 



df) 



—k 



B 



'l + a)cT a V In 



v 

cT a V 
n + cT a V 



+ n In 



71 



Helmoltz free energy F 



F = k B T 
Internal energy U 



n In 



71 



n + cT a V 



+ cT a V In 



n + cT a V 



cT a V 
n + cT a V 



(25) 



(26) 



U = —ksT 2 ln(l - z) 



-k B T ak ln(l - z) = a k B TN; (27) 



Specific heat at constant volume Cy 

'n + cT a V 



C v = k B a(a + l)cT a V In 



cT a V 



71 



'n + cT a V 



a 2 cT a V. (28) 



By comparing now the internal energy of an ideal gas U = f(u) k B T In Q 
(/(f) is here a function of the degrees of freedom of the system under investi- 
gation) with the result of our calculation (27) and by recalling that N = In Q 
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it follows that a = f(v), i.e. a is a function of the clan degrees of freedom v. 
Notice that clan aggregation process is determined by the increase of the 
clan degrees of freedom, i.e. as the collision becomes harder clan degrees of 
freedom should be larger and a phase transition from a class of events to 
another one can occur by varying clan degrees of freedom. 



5 A phase transition between the semihard 
and hard component at parton level implies 
a discontinuity in the average number of 
particles of the two components at hadron 
level 

Within the framework discussed in the previous Section we would like to 
show: 

I. that the expected phase transition in the ideal clan gas is a first order 
phase transition and that it is revealed at hadron level through a sudden dis- 
continuity in the average number of final charged particles (the counterpart 
of the remark given in Ref . [8] that Pn vs. n plot for the semihard class of 
events from concave k semihard > 1 is becoming convex with khard < 1); 

II. that in the maps of zeroes of the truncated NB (Pascal) MD grand canon- 
ical partition function in the rescaled complex w-plane u = 1 is an accumu- 
lation point of zeroes (from inside for k semi h a rd > 1 and from outside for 

hard < 1 [ZD- 
Assume then that the partonic clan gas of the semihard component in pp col- 
lisions has below the critical cm. energy, ^/s~ C) ^semihard degrees of freedom 
and above yfsl the partonic clan gas of the hard component has Uhard degrees 
of freedom, with of course a hard = f{v hard ) > a semihard = f{v 'semihard)- Re- 
calling now that the truncated grand canonical partition function is analytic, 
the phase transition from the semihard to the hard ideal clan gas is expected 
to occur at parton level (p) in the thermodynamic limit, corresponding to 

— (semihard) — (hard) 

l im ™P = (semihard) = ^ % _ (hard) 

(semihard) rAsemihard) re _( ha rd) Tr(hard) "c 

n p '^oo Vp rip '— too Vp 

T ,( semihard) ,.{hard) 

V~i ; ^co V„ '— too 
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i (semihard) -i (hard) ... i -i ., . £ 

where p c and p c are critical parton densities of the two phases. 

In addition being the Helmoltz free energy density for the two phases 

f(0= i im 

with i = semihard or /iard, at the critical temperature T c one has 

r (semihard) r(hard) 

J p J p 



i.e. 



which leads to 



semihard __ ^ ^^hard 
^■semihard C ^hard C 



Qf,, , \ a hard a semihard 

^semihard 1 



C OL hard 



rj~> J ^semihard J ^29^ 



Since c a .. , and c Q , , do not depend on temperature T and volume V, 

^semihard <-*-nard -L ' 

they are expected to be the same at hadron (h) and parton (p) level and by 
recalling Eq.(24), it follows that 



T c = lim 



rp°> semihard y(semihard) 
h h 



^hard a semihard 



, .{semihard) 

V h \ rp a hardy(hard) 



k 



yihard)^ 



Tr(hard) \ 



*hard a semihard 



^ h , A ( se miWd) I ^(semihard) J 

A relation which by assuming that the particle densities are continuous also 
at hadron level, 

— (semihard) — (hard) 

l im U h = (semihard) = ^ Uf, = ( hard) 

(semihard) r(semihard) r H (hard) -,r(hard) ' n 

n h V h n h V h 



, A semihard) , J hard) 



implies 



i 



— (hard) \ a hard a semihard 



_ (semihard) 

_(hard) 
h 



Tc = Th lim | -p-^yj . (30) 



It follows from T c > T ft that Eq.(30) is certainly satisfied when 



— (hard) 
"h 



n 



!»«.. 7(=w) > 1- (31) 



_ {semihard) 



_ (hard) 
h 
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Therefore the phase transition from the semihard ideal clan gas to the hard 
one, occurring at parton level, is revealed at hadron level by a discontinuity 
in the average number of final particles of the two clan gas components. 
By using GLPHD Eq.(31) can of course be rewritten as follows 

-(hard) -.(hard) 

, A. > m 

_ (hard) 

rip -^>oo 

but since at the phase transition one has 

-(semihard) -(hard) 

U p — IL p 



we get 

^.(hard) 



> 1 (33) 



^■(semihard) 

i.e. 

^.(hard) ^ ^.(semihard) (34) 

Therefore Eq.(30) becomes 

( (hard) \ a hard -a 3emihard 
t — "7T (35) 
^■(semihard) I y ' 

suggesting that hadronization process differs in the two classes of events, 
under the assumption that the hadronization temperature is the same. 
In addition, being in the thermodynamic limit the entropy density y equal 
to 



S ^ 
V = -k B 



(1 + a) cT a In I — < ^——- \ + p c In ^ f> ' 



Pc + cT a J ^ \ Pc + cT a 



(36) 



it is clear that 

lim — ^ lim — 

i.e. the phase transition from the semihard to the hard partonic ideal clan 
gas corresponds (being S — — (j/?) nV ) to a first order phase transition. 

Notice that the same conclusion can be achieved by studying the zero 
distribution properties of the M-truncated grand canonical partition function 

Q M = l + zZ x + z 2 Z 2 + ... + z M Z M . (37) 



16 



From a known algebraic relation between coefficients and zeroes of a polyno- 
mial (see Ref.[15j). in the thermodynamic limit (M, V —>■ oo with M/V = p) 
at the critical temperature T c one gets 



M 

n=l 

M 

n=l 



J M-\ 



M M M,V^oo p c 



Z M M + k-1 M + cT a V-l Pc + cT< 



M / ry \ 2 ry 



Zm J Z M 

' M\ 2 



M,V^oo ( p. 



(f + cT»-i) (f + cT»-i) 2 Vp. + cT, 



M / \ M 

M M,V->oo / p, 



n=l 



E4' -(^) (38) 



In Appendix it is shown that all the zeroes of a degree M polynomial are 
distinct and lie on the circle of I/7 radius at angles 6 n = n^^j when 



- / 1 V 

J2 z n = -(-J > J = M,...,M. 



Notice that in the present approach, I/7 is given by 

I Pc 

7 ~ p c + cT«' 



(39) 



which corresponds to the critical fugacity z c and to a zeroes accumulation 
point, i.e. to a non-analicity point of the grand canonical partition function 
(4). Now, since the zeroes approach the positive real axis with a 7r/2 angle, 
it is confirmed that the phase transition is of first order [H]. In fact it can be 
shown that the linear density of zeroes on the positive real axis is different 
from zero and equal to 

1/M , 1 M + 1 1 
lim — — = lim 



M-foo z c 9\ m-*oo M 2nz c 2nz c 
Eq.(29) is obtained by using the identity 

Pc Pc 



n _|_ „ rjiOi semihard n I „ H^^hard ' 

Pc T <-« 3eml hard J- c p c T ^ct hard -Lc 
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i.e. by asking that I/7 (the zeroes accumulation point) is the same for the 
two phases. 

Notice that the properties of zeroes of the NB grand canonical partition 
function in the thermodynamic limit, Eq.s(38), do not depend on the shape of 
temperature dependence of parameter k. Therefore, the linear dependence of 
parameter k on volume V, as described by Eq.(23), contains two information: 
1) the presence of a phase transition (the zeroes stay on the circle of radius 
equal to the critical fugacity z c , which corresponds to a zeroes accumulation 
point) and 2) the phase transition is of first order (since the zeroes stay on a 
circle, they approach the positive real axis with a 7r/2 angle |HJ). 
Notice that the general result given by Eq.(39) in the limit p c —>■ 00 implies 
z c > 1 , as shown in Ref. [7j. 

Appendix: On a theorem on the zeroes of a 
degree M polynomial 

Let us consider a generic polynomial of degree M: 

M 

G M (z) = zmp m = P o + zPi + z 2 ?2 + ... + z M P M . (40) 

m=0 

We now prove the following theorem: if the zeroes of the algebraic equation 

G M {z) = (41) 

satisfy the conditions 

M /<y 



M / 1 \ j 

$>n = - - ' 3 = 1 > 2 >~> M (42) 



n=l 



2tt 



they are all distinct and lie on the circle of I/7 radius at angles 9 n = n M , j . 
with n = 1,2, ...,M. 

Proof: since the zeroes of a polynomial occur in complex conjugate pairs, 
Eq.s(42) can be rewritten as follows 

M M 



n=l 

- f 1 \ j 

J2\ z n\ J cos[je n ] = -(-j , j = l,2,...,M. (43) 



n=l n=l 

M 
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The request that the zeroes of Eq.(41) lie on the circle of radius I/7, implies 
of course that | z n \— I/7, Vn = 1,2, ...,M, and therefore the following 
condition for the 9 n angles holds: 

M 

]Tcos[j^] = -l, j = l,2,...,M. (44) 



n=l 



The theorem can be proved simply by showing that Eq.s(44) are indeed 
verified, as will be shown in the following. Notice first that 



M 



$^cos[> 



nx\ 



n=0 



3? 



ft 



M 



3? 



J(M+l)x _ 
e ix _ I 



n=0 

e i(M+l)x/2 _ e -i(M+l)x/2 gi(M+l)rr/2 



D ix/2 



-ix/2 



D ix/2 



sin[(M + l)x/2] ^ [e!Ml/2] = cos[Mx/2] sin[(M + l)x/2] 



sin[x/2] 



sin[x/2] 



and of course 



M 

E 

n=l 



cos 



cos[Mx/2] sin[(M + l)x/2] 
sin[x/2] 



Calling now nx = n , with 6 n = nj^, one gets 



M 



M 



^cos[# n ] = 



COS 



71=1 



2tt 



n- 



M+l 

n=l L 

COS [ M m] Sill [ 7r ] 

sin [l^T] 



oosgM^] sin[|(M + 1) 



2w 1 
M+lJ 



1 = 



1 = -1. 



Accordingly, by writing ra as with n = n^jr^ and j = 1,2, ...,M, it 
follows that 

M 

^cos[j6> n ] 



M 

E 



cos 



n=l 



2tt 



M+ 1 



cos[|M^]sin[|(M+l) 



2tt 1 
M+l-l 



n=l 

cosfj'M-^] sin[j7r 



sin[i^ 



- 1 



2 M+lJ 



- 1 = -1, 



and therefore 



M+lJ 



E 

n=l 



COS 



2tt 



M+ 1 



= -1, 



with j = 1,2,...,M. Since the set of zeroes of a polynomial is univocally 
defined, the theorem is proved. 
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